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We calculate cohesive properties of LiF, NaF, KF, LiCl, 
NaCl, and KCl with ab-initio quantum chemical methods. 
The coupled-cluster approach is used to correct the Hartree- 
Fock crystal results for correlations and to systematically im- 
prove cohesive energies, lattice constants and bulk moduli. 
After inclusion of correlations, we recover 95-98% of the total 
cohesive energies. The lattice constants deviate from exper- 
iment by at most 1.1%, bulk moduli by at most 8%. We 
also find good agreement for spectroscopic properties of the 
corresponding diatomic molecules. 



I. INTRODUCTION 

One of the earliest methods for a quantitative treat- 
ment of the cohesion of ionic solids was the Born-Mayer 
theoryEI. Lowdinfl made a first quantum-mechanical 
approach starting from the symmetrically orthogonal- 
ized orbitals of the free ions; these orbitals were used 
to approximate the density matrix and to calculate 
the Hartree-Fock energy. Since the advent of density- 
functional theory and especially the local density ap- 
proximation, the latter methods have become standards 
in solid state physicaj. However, there has also been 
progress in the development of wavefunction based meth- 
ods. Hartree-Fock (HF) calculations can be done rou- 
tinely nowadays with the help of the program package 
Crystalo, and it is even possible to include electron cor- 
relations. One way of achieving that is by multiplying the 
HF wavefunction with a Jastrow factor containing several 
parameters; these parameters, can be optimized with the 
help of Monte-Carlo methodsO. A first attempt to include 
correlations by means of quaiiUim chemical methods was 
made using the Local AnsatzQu; here local excitation op- 
erators are applied for modifying the HF wavefunction. 
In the last years, an "incremental scheme" (an expansion 
of the total correlation energy in terms of one-body, two- 
body, three-body and higher contributions, the so called 
"local increments") has been developed and successfully 
applied to semiconductoraa. This method has been ex- 
tended to iwiic solids and applied to several oxides (MgO, 
CaO, NiO)Ll. Alkali halides are model examples of ionic 
solids and haye recently been carefully investigated at 
the HF levelllJ. The major part of the experimental lat- 
tice energy is already recovered at this level. However, 



the lattice constants significantly deviate from the experi- 
mental values, especially for the heavier compounds. We 
want to show that the incremental scheme can explain 
the deviations of the HF results from experiment. 



II. THE METHOD 

A. Incremental scheme 

The scheme has been explained in earlier work and we 
only repeat the main ideas. The correlation energy of 
the solid is expanded into a sum of local contributions 
(increments) 



A A,B 

+ 1 ^ Ae{ABC) + ... 



A,B,C 

e{A) is the correlation energy of a group of local- 
ized orbitals (a so-called one-body increment), the non- 
additivity Ae{AB) = e{AB) - t{A) - e{B) defining a 
two-body increment, and so on. Usually, this series is 
evaluated up to three-body increments. The increments 
are extracted from clusters containing up to three explic- 
itly described ions (i.e. ions with a high-quality basis set) 
embedded in a set of pseudopotentials and point charges. 
They should be well transferable, which means that they 
should only weakly depend on the specific cluster cho- 
sen for their evaluation (e.g. the value of a one-body 
increment obtained from a cluster with one explicitly 
described ion only weakly varies when extracted from 
a cluster with more than one explicitly described ion). 
As correlation scheme, we chose the coupled-cluster ap-, 
proach with single and double substitutions (CCSD)tll 
with an exponential ansatz for the correlated wavefunc- 
tion: 



/ 



|*CCSD >= exp 



Ers + + 



|*SCF > 



In addition, we appHed the CCSD(T)il scheme includ- 
ing triple excitations in a perturbative way. All the cal- 
culationSpjKace done with the ab-initio program package 
MOLPRCOO. LocaHzation was done by the Foster-Boys 
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methodEi, and all of the ns,np valence and outer-core 
orbitals of the halide and alkali ions, respectively, (n = 
2 for F, Na and 3 for CI, K, Is in the case of Li), were 
correlated. 



B. Pseudopotentials and basis sets 

The increments are taken from cluster calculations. 
The ions to be correlated are accurately described with 
extended basis sets. Negatively charged ions are embed- 
ded with X+ pseudopotentials as next neighbors to simu- 
late the Pauli repulsion. Finally, the system is embedded 
in a set of point charges (typically 7x7x7 lattice sites 
with charges ±1 in the interior and reduced by factors 
of 2, 4 and|-8. at the surface planes, edges and corners, 
respectivelyllj) . The description of the expliciily treated 
ions is as follows. We used a [5s4p3(i2/] basiaiZI for F and 
a [6spfl3d2/] basifO for CI. For Li, we used a [5s4p3d2/] 
basiill, for Na a [7s6p5d4/] basis (Ref. with d and 
/ functions uncontracted) . Finally, for K we used a 9- 
valence-electron pseudopotentiaLJ with the correspond- 
ing sp basis set (^contracted) and augmented with 5 d 
and 3 / functionsEj, resulting in a [7s6p5c?3/] basis. 



III. RESULTS 

A. Ionization potentials, electron affinities and 
results for the diatomic molecules 

In Table |, we give results for atomic electron affini- 
ties and ionization potentials. At the correlated level, we 
obtain good agreement with experiment (to <0.1 eV) in 
all cases. Results for the diatomic molecules are given 
in Table Q. Again, we obtain nice agreement, to < 
0.02 A (1%) for bond lengths, 24 cm^i (4%) for vi- 
brational frequencies, and 0.1 eV for dissociation ener- 
gies De- Note that we calculated Df. as the difference 
Eatomi+ Eatom2- Ediatomic, in contrast to Rcf. |2l|, where 
the dissociation energy was first calculated with respect 
to the singly charged ions and then corrected with the 
help of the experimental electron affinities and ionization 
potentials. The experimental dissociation energy for NaF 
from Ref. |2^ isprobably too high, the experimental value 
given in Ref. ^ (Z3e=4.97 eV) and the theoretical value 
from Ref. |2ll are closer to our calculated value. 



B. Results for the solid 

1. Hartree-Fock calculations 

We repeated the Crystal calculatipns from Ref. |l^ 
with essentially the same basis setso. We calculated 
both the lattice energy (cohesive energy with respect to 
the ions) as well as the cohesive energy with respect to 



the neutral atoms. The lattice energy is already in good 
agreement with experiment. This is what one would ex- 
pect since in purely ionic solids (the MuUiken population 
analysis gives a charge transfer very close to ±1 in all 
cases) the Madelung energy makes the most important 
contribution to the lattice energy; the Madeluag energy is 
already in rough agreement with experimentfij. However, 
the cohesive energy with respect to the atoms is less well 
described as a consequence of the missing intra-atomic 
correlation effects. Moreover, lattice constants are by up 
to ^ 5% too large at the HF level, bulk moduli up to ~ 
21% too small. 



2. One-body increments 

Results for the crystal correlation energies are given 
in Tables II],rv|,[v|,[V^. Concerning the one-body incre- 
ments, we obtain nearly the same correlation energy for 
the free alkali ions and the corresponding embedded ions. 
This is of course a consequence of the small ionic radii 
of the cations. In the case of the anions F^ and Cl^, 
we find that the absolute value of the correlation energy 
in the solid is smaller than for the free ion, by up to 0.4 
eV. Such an effect was already found in the calculations 
on the oxideaJ and is explained by the lower level spac- 
ing of the excited states for the free ion compared to the 
embedded ion where excitations are higher in energy. 



3. Two-body and three-body increments 

The two-body correlation-energy increments decrease 
rapidly. The decay is compatible with a van der Waals 
law from second nearest neighbors on, cf. Table By 
far the largest contributions come from next-neighbor 
metal-halide (M-X) and halide-halide (X-X) interactions. 
The total effect of the M-X inter-atomic correlations is 
similar for X = F and X = CI, but for given X increases 
from Li to K (i.e. with increasing polarizability a of the 
metal ion) in such a way, that the ratio of the M-X con- 
tribution to the X-X contribution changes from < 1 to 
> 1 (cf. Tables 1^, 0). The X-X increments in turn are 
larger in magnitude for CI than for F, in agreement with 
the trend of the respective a values but in contrast to the 
situation for the intra-atomic difference in correlation en- 
ergies e(free ion) — e(embedded ion). Quantitatively com- 
paring the F-F and Cl-C l nex t-neighbor increments from 
different systems (Table VII ) and assuming a purely van 
der Waals interaction, we find that even in that case the 
van der Waals-law holds surprisingly well. The Cq coef- 
ficient can be determined from the two-body increments. 
For the sake of simplicity, we assume a purely van der 
Waals interaction already for next neighbors and for all 
types of coixelations (e.g. also spin-flip processes for Ni-0 
incrementscl) . The result for Cq — AE x obtained this 
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way is comparable to results from literature, e.g. Refs. 
^ and references therein. 
An estimate of the van der Waals interaction can 
be obtained^using the London formula for dispersion 
interactionscS: E = ~ /p^^+fp^ with the ioniza- 
tion potentials (IP) as characteristic excitation energies 
and polarizabilities (a) of the two interacting systems 
(ji is of order unity, r is the distance). Polarizabili- 
ties and ionization potentials were calculated with the 
same arrangement as the one-body increments: One ion 
with extended basis set was embedded in a set of point 
charges at the experimental lattice constant (and pseu- 
dopotentials as next neighbors, in the case of anions) . To 
evaluate the polarizabilities, we applied a small dipolar 
field and find .values in good agreement with values from 
literatureCJ'HEd. The ionization potential was calculated 
with the same cluster, which is certainly a crude approx- 
imation because effects such as long-range polarization 
are not included: the IP obtained this way is not what 
would be experimentally measured for the solid. Our 
CCSD results for the two-bo dy i ncrements are roughly 2 
to 5 times larger (see Table VII) than what is predicted 
from the London formula. This implies that the London 
formula can give a qualitative understanding of the mag- 
nitude of the interionic interaction and the parameters 
describing it {a, excitation energies), but is not able to 
predict results quantitatively. Van der Waals interactioBS. 
in extended systems hav£,also been considered for HetHi 
(see also a recent reviewE3). 

We calculated three-body increments only for KCl (Ta- 
bles and 0). We find that they are very small indi- 
cating a rapid convergence of the incremental expansion. 
Neglecting three-body increments is not a serious approx- 
imation, therefore. 



4- Sum of increments and discussion 



The sums of the increments are given in Table III . In- 
cluding correlations, we obtain 95 to 98% of the experi- 
mental cohesive energies. The relatively good agreement 
of the HF lattice energies already mentioned above turns 
out to be due to a partial error cancellation. When the 
HF cohesive energies are calculated with respect to the 
free ions, the corrections due to the missing correlation 
effects have opposite signs: the one-body contributions 
diminish the cohesive energy since the absolute value for 
the free anion is higher than that for the embedded ion; 
on the other hand, the van der Waals interactions which 
are also missing at the HF level lead to an increase of 
the cohesive energy (cf. Tables 0, 0). The compen- 
sation is nearly perfect for LiF, but already for KF the 
inter-atomic correlation effects overcompensate the intra- 
atomic ones by nearly a factor of 2, and the weight is still 
further shifted in favour of the two-body effects for the 
MCI crystals, so that for KCl, e.g., a factor of ^6 is 
reached. 



After inclusion of correlations, the lattice constants de- 
viate by at most 1.1% from experiment. As already found 
in the context of the oxides, the one-body increments 
would enforce larger lattice constants (the absolute value 
of the correlation energy of an anion increases when the 
lattice constant increases because of the lower level spac- 
ing at larger lattice constant) . The large reduction of the 
lattice constants, on the other hand, is a two-body effect 
resulting from the van der Waals interaction between the 
ions. The CC|Sp(T) results turn out to be slightly supe- 
rior to ccsdeI. 

At fixed lattice constant, inclusion of correlations leads 
to a decrease of the bulk modulus. However, for most of 
the solids considered here correlations reduce the lattice 
constant. This means that the HF bulk modulus has to 
be calculated at a smaller lattice constant where it in- 
creases again. As a net result, correlations increase the 
bulk moduli in most cases. Note that the bulk moduli 
are more sensitive to the fitting procedure than cohesive 
energies and lattice constants and that they also have 
large experimental uncertainties even at room tempera- 
ture (see the comparison in Ref. 34). 

A more detailed account of correlation contributions 
to the potential-energy surface of KCl is given in Fig. 
|l|, where we display the difference of correlation energies 
e(embedded Cl~) -e(free Cl~) as a function of the lattice 
constant, i.e. its variation from free Cl^ to an embedded 
Cl^ in KCl. Starting from a very small (unrealistic) lat- 
tice constant a, the correlation energy e(embedded CI ) 
decreases in magnitude with increasing a - excitations 
into dxy, dyz, c?a;z-orbitals are very important for small a 
since these orbitals have smaller overlap with the region 
that is occupied by the K electrons -, then passes through 
a minimum and increases again because of the argument 
given earlier (excitations into the diffuse CI 4p orbitals 
are lower in energy the larger the distance to the K elec- 
trons). The next-neighbor K-Cl and Cl-Cl correlation- 
energy increments also shown in Fig. |l| monotonously 
decrease with increasing a, for larger distance according 
to the van der Waals law. The three contributions de- 
picted in Fig. ^ are the most important ones and nearly 
exhaust the incremental expansion (see Table VI , the re- 
maining increments amount to ~ 1 mH only). The first 
derivative of their sum with respect to the lattice con- 
stant clearly shows that in total correlations reduce the 
lattice constant. The second derivative shows that - at 
fixed lattice constant - correlations reduce the bulk mod- 
ulus (the one-body increments alone might lead to an 
increase of lattice constant and bulk modulus, but are 
outweighed by the two-body increments). 

Several density functional calculations are avail- 
able in literatum. for the systems considered. A 
KKR calculationEEl |j(combined with a local exchange- 
correlation potentiaJcZI) and more recently a full-potential 
xc-LDA calculationEj have been performed. In Refs. |3|, 
I39I correlation-only density functionals with gradient cor- 
rections have been included a posteriori (i.e. using the 
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density and non-local exchange energy from a Hartree- 
Fock calculation). The best density functional results 
are in good agreement with experiment, but it seems to 
be difficult to choose one single functional as reference 
method. 

In Rcf. a large number of alkali halide clusters has 
been investigated. Bulk properties were extrapolated 
from cluster calculations by linearly fitting the energy 
vs. n~^/'^, where n is the number of MX units. The re- 
sults for the lattice energies Eiat are in good agreement 
with experiment. The predicted correlation corrections 
are in agreement with our findings for LiF (~ 0), but 
different for NaCl (an increase of \Eiat \ of ^ 0.003 H is 
reported, we find ~ 0.013 H at the CCSD level) and KCl 
{A\Eiat\ ~ -0.011 H from Ref. we obtain a CCSD 
value of ~ 0.016 H). The geometries were optimized at 
the HF level using a M32N32 cluster. It was proposed to 
use the bond length of the interior cube of this cluster 
as an estimate of the lattice constant of the solid. This 
leads to a slight underestimation in all cases compared 
to the HF lattice constants from Crystal calculations. 
Surface effects are probably the explanation for the dif- 
ferences, since each atom of the interior cube has three 
next neighbors also residing in the interior cube, but also 
three next neighbors located at the surface whose charges 
will be different from interior ions; the Pauli repulsion 
and the Madclung field arc probably not too well repro- 
duced. This is avoided in our approach since a cluster 
approach is applied at the correlated level only, and even 
there all explicitly treated ions are surrounded by pseu- 
dopotentials (or point charges) simulating bulk cations 
(or anions). 

IV. CONCLUSION 

We have shown that the method of local increments 
can successfully be applied for the determination of bulk 
electron-correlation effects in alkali halides. The main 
shortcoming of the Hartree-Fock approximation is the 
missing inter-ionic van der Waals interaction which re- 
sults in too large lattice constants (by up to 5%). Af- 
ter including correlations at the coupled-cluster level, the 
deviations of the lattice constant from the experimental 
values are reduced to a maximum of 1.1%. We obtain be- 
tween 95 and 98 % of the cohesive energies with respect 
to neutral atoms or 97 to 98 % of the lattice energies. 
Bulk moduli exhibit satisfactory agreement with experi- 
ment, with a maximum deviation of ~ 8%. 
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FIG. 1. CCSD correlation-energies for KCl as a function 
of the lattice constant. The two-body increments are al- 
ready multiplied with the corresponding weight factors. Dis- 
played is the difference of correlation energy e(embedded Cl~) 

-e(free Cl~) (dashed ), the two-body increment Cl-Cl for 

next neighbors (dashed-dotted - • - ), the two-body increment 
K-CI for next neighbors (dotted • ■ •), and the sum of these 

three correlation energies (solid line ) which make the 

most part of the total correlation contribution to the cohesive 
energy. 
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TABLE 1. Electron affinities and ionization potentials in 
Hartrees. (1 H = 27.2114 eV) 



System 


HF 


CCSD 


CCSD(T) 


exp 


11 


F ^ F" 


0.05070 


0.11612 


0.12192 


0.12499 


ci cr 


0.09505 


0.12605 


0.12919 


0.13276 


Li ^ Li+ 


0.19631 


0.19731 


0.19733 


0.19814 


Na ^ Na+ 


0.18195 


0.18785 


0.18810 


0.18886 


K ^ K+ 


0.14679 


0.15637 


0.15723 


0.15952 
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TABLE II. Bond lengths Re (A), dissociation energies 
De (eV) and vibrational frequencies Ue (cm"'^) of diatomic 
molecules. The values taken from literature are CI(SD) cal- 
culations. 





RHF 


CCSD 


CCSDI'T'I 


literature 


21 


expt 


22 


LiF 














Re 


1.555 


1.561 


1.565 


1.571 


1.564 




943 


923 


910 


919 


910 


D, 


4.12 


5.85 


5.98 


6.12 


5.97 


NaT 














Re 


1.924 


1.925 


1.929 


1.921 


1.926 




549 


517 


512 


538 


536 




3.11 


4.77 


4.91 


5.02 


5.36 


KF 














Re 


2.204 


2.189 


2.189 


2.184 


2.171 


LJe 


420 


422 


421 


428 


428 


De 


3.29 


4.89 


5.03 


5.10 


5.10 


LiCl 














Re 


2.037 


2.026 


2.028 


2.033 


2.021 


LJe 


645 


645 


642 


646 


643 


De 


3.85 


4.76 


4.84 


4.86 


4.8t 


i 


NaCl 














Re 


2.390 


2.344 


2.344 


2.366 


2.361 


UJe 


359 


368 


367 


361 


366 


De 


3.26 


4.12 


4.20 


4.23 


4.25 


KCl 














Re 


2.738 


2.692 


2.688 


2.697 


2.667 


Ue 


266 


276 


276 


273 


281 


De 


3.48 


4.21 


4.29 


4.33 


4.36 
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TABLE III. Hartree-Fock (HF) and correlated results (CCSD, CCSD(T)), in comparison to density-functional (DFT) and 
experimental values, for the solids. Cohesive energies E (with respect to neutral atoms) and lattice energies Eiat (with respect 
to free ions) are given in Hartree units, lattice constants a in A and bulk moduli B in GPa. Zero point energies have been 
estimated with a Debye approximation (Debye temperatures taken from Ref. |4^ ) and added to the experimental cohesive 
energies. The experimental bulk moduli are at 4.2 K and have been taken from Ref. k3 and references therein. 





HF 


CCSD 


CCSD(T) 


DFT 


exptCil'Htill 


LiF 












Elat 


0.3975 


0.3976 


0.3961 


0.417", 0.400^ 0.365^= 


0.404 


E 


0.2534 


0.3179 


0.3222 


0.352'', 0.345 " 


0.331 


a 


4.011 


3.991 


3.993 


4.035 4.05 " 3.88", 3.96\ 4.13"^ 


4.010 


B 


78.9 


70.1 


74.9 


78.3 70.5 ^ 95 ", 83 \ 60 " 


69.9 


NaF 












Elat 


0.3496 


0.3518 


0.3504 




0.358 


E 


0.2186 


0.2803 


0.2845 


0.323 0.294 " 


0.294 


a 


4.636 


4.601 


4.603 


4.582 4.76 " 


4.609 


B 


52.2 


55.7 


53.9 


55.8 42.3 " 


51.4 


KF 












Elat 


0.3028 


0.3101 


0.3100 




0.318 


E 


0.2076 


0.2707 


0.2755 


0.294 " 


0.283 


a 


5.450 


5.331 


5.320 


5.40 " 


5.311 


B 


29.9 


34.4 


34.8 


31.3 " 


34.2 


LiCl 












Elat 


0.3088 


0.3225 


0.3241 




0.331 


E 


0.2096 


0.2533 


0.2580 


0.251 0.265 " 


0.266 


a 


5.281 


5.136 


5.124 


5.32 5.08 " 


5.106 


B 


30.1 


35.2 


34.8 


28 35.2 " 


35.4 


NaCl 












Elat 


0.2839 


0.2960 


0.2971 


0.304", 0.312\ 0.285=, 0.307-'', 0.303^, 0.300'''' 


0.302 


E 


0.1978 


0.2350 


0.2390 


0.239^*, 0.232*'' 


0.246 


a 


5.791 


5.646 


5.634 


5.737 5.75 ^ 5.47", 5.49\ 5.83^ 5.53^, 5.51^ 5.54'^-^ 


5.595 


B 


24.5 


26.6 


26.6 


25.5 22.8 31", 29\ 21^ 32.5^, 32. 1», 30.1'*'' 


26.6 


KCl 












Elat 


0.2538 


0.2687 


0.2704 




0.275 


E 


0.2035 


0.2398 


0.2438 


0.249 0.243 " 


0.248 


a 


6.548 


6.314 


6.295 


6.30 6.26 " 


6.248 


B 


15.5 


18.4 


21.3 


19.7 18.9 


19.7 



°Rcf. 38 
^Ref. 38 
'^Ref. 38 
36 



34 

i 

39 
38 



Hartree-Fock exchange, Perdew and Wang 91 correlation functionalcl 
LDA exchange and correlation l— , 
Becke exchangee^ and Perdew and Wang 91 correlatip« functionala 
KKR calculation with local exchange and correlationcZl 
LDA exchange and correlation ■— ■ 
Hartree-Fock exchange, Colle and Salvetti correlation functionalEa 
Hartree-Fock exchange, Perdew 1986 correlation functionalcZl ■— ■ 
Hartree-Fock exchange, Perdew and Wang 91 correlation functional^ 



*For further density functional results for NaCl, see also Ref. 
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TABLE IV. Local correlation-energies per primitive unit cell (in Hartree) for LiF ( at a lattice constant of 3.99 A), NaF 
(4.60 A), and KF (5.34 A) 





LiF 

CCSD 


CCSD(T) 


NaF 

CCSD 


CCSD(T) 


KF 

CCSD 


CCSD(T) 


free X+ ^ embedded X+ 


-0.000021 


-0.000021 


-0.000165" 


-0.000179" 


-0.000013 


-0.000016 


free F~ — > embedded F~ 


+0.011776 


-h0.014782 


+0.010106 


+0.012820 


+0.009684 


+0.012352 


sum of F-F increments 


-0.007926 


-0.009141 


-0.003384 


-0.003954 


-0.001170 


-0.001359 


sum of X-F increments 


-0.003970 


-0.004254 


-0.008554 


-0.009346 


-0.014940 


-0.017074 


sum of X-X increments 


-0.000018 


-0.000018 


-0.000198 


-0.000210 


-0.001422 


-0.001605 


sum 


-0.000159 


+0.001348 


-0.002195 


-0.000869 


-0.007861 


-0.007702 



° See footnote (Ref. ||). 

TABLE V. Local correlation-energies per primitive unit cell (in Hartree) for LiCl (at a lattice constant of 5.14 A), NaCl 
(5.65 A), and KCl (6.30 A). 





LiCl 






NaCl 


KCl 






CCSD 


CCSD(T) 


CCSD 


CCSD(T) 


CCSD 


CCSD(T) 


free X+ embedded X^ 


-0.000013 


-0.000013 


-o.oooior 


-0.000109" 


-0.000005 


-0.000005 


free Cl~ embedded CI" 


+0.002567 


+0.003572 


+0.002448 


+0.003415 


+0.002426 


+0.003411 


sum of Cl-Cl increments 


-0.014439 


-0.016785 


-0.008124 


-0.009495 


-0.003732 


-0.004368 


sum of X-Cl increments 


-0.002712 


-0.002906 


-0.007112 


-0.007746 


-0.014992 


-0.017066 


sum of X-X increments 


absolute value <10 ^ 


-0.000054 


-0.000060 


-0.000444 


-0.000501 


sum of three-body increments 










+0.000388 


+0.000372 


sum 


-0.014597 


-0.016132 


-0.012943 


-0.013995 


-0.016359 


-0.018157 



" See footnote (Ref. ||). 
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TABLE VI. Local correlation-energies per primitive unit 
cell (in Hartree) for KCl at a lattice constant of 6.57 A. The 
quantities involving 2 and 3 ions are non-additivity correc- 
tions (increments). 





weight 


('('SD 


CCSD(T) 1 


free embedded 


1 


-0.000004 


-0.000004 


free CP ^ embedded CP 


1 


+0.002059 


+0.002911 


C1(0,0,0)-C1(0,L1) 


6 


-0.002736 


-0.003228 


C1(0,0,0)-C1(2,0,0) 


3 


-0.000138 


r\ r\r\r\-\ 

-0.000162 


C1(0,0,0)-C1(2,1,1) 


12 


-0.000144 


-0.000168 


C1(0,0,0)-C1(2,2,0) 


6 


-0.000030 


-0.000036 


K(0,0,0)-C1(1,0,0) 


6 


-0.011256 


-0.012858 


K(0,0,0)-C1(1,1,1) 


8 


-0.000320 


-0.000360 


K(0,0,0)-C1(2,1,0) 


24 


-0.000192 


-0.000216 


K(0,0,0)-K(0,1,1) 


6 


-0.000318 


-0.000360 


K(0,0,0)-K(2,0,0) 


3 


-0.000018 


-0.000021 


C1(1,0,0)-C1(0,1,0)-C1(0,0,1) 


8 


-F0.000064 


+0.000080 


C1(0,0,0)-K(0,1,0)-C1(0,1,1) 


12 


-1-0.000204 


+0.000204 


sum 




-0.012829 


-0.014218 
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TABLE VII. Comparison of CCSD two-body increments AE between next neighbors (without multiplying with the weight 
factor). All results are given in atomic units (except for the lattice constant in column 2). r is the distance between the 
respective ions in bohr. 



System 


lattice 
constant a in A 


AE 


AE X 


IPcat 


IP 

an 


Ocat 




2 Vi'^ X AE 

U <^IL\<-IL'2 11 L 1 ^ 2 


F-F (LiF) 


3.99 


-0.001181 


-27.1 


2.3 


0.52 


0.19 


5.0 


2.8 


F-F (NaF) 


4.60 


-0.000502 


-27.1 


1.3 


0.47 


0.97 


5.4 


2.6 


F-F (KF) 


5.34 


-0.000174 


-23.0 


0.80 


0.42 


5.4 


5.4 


2.5 


Cl-Cl (LiCl) 


5.14 


-0.002155 


-226 


2.4 


0.45 


0.19 


19 


1.9 


Cl-Cl (NaCl 


5.65 


-0.001215 


-225 


1.4 


0.42 


0.97 


19 


2.0 


Cl-Cl (KCl) 


6.30 


-0.000558 


-199 


0.85 


0.39 


5.4 


18 


2.1 


O-O (MgO)" 


4.18 


-0.002582 


-78.4 


2.1 


0.38 


0.48 


9.7 


2.9 


O-O (CaO)" 


4.81 


-0.001067 


-75.2 


1.2 


0.27 


3.1 


9.7 


3.9 


O-O (NiO)" 


4.17 


-0.003356 


-100 


0.42 


0.41 


2.8 


11.4 


2.5 


Li-F 


3.99 


-0.000627 


-1.80 










3.0 


Na-F 


4.60 


-0.001351 


-9.11 










3.4 


K-F 


5.34 


-0.002382 


-39.3 










3.3 


Li-Cl 


5.14 


-0.000440 


-5.77 










2.8 


Na-Cl 


5.65 


-0.001132 


-26.2 










2.9 


K-Cl 


6.30 


-0.002392 


-106 










2.7 


Mg-O" 


4.18 


-0.003129 


-11.9 










5.3 


Ca-0" 


4.81 


-0.005906 


-52.0 










5.2 




4.17 


-0.009958 


-37.3 










3.8 



Ref. I 
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